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HARMONIC  SYNTHESIZER 
HAVING  20  SINE  AND  20  COSINE  TERMS 

During  work  on  structural  transformations  in  crystals  we  have 
had  the  problem  of  determining  the  structures  of  the  phases  Involved. 
In  this  connection,  determination  of  the  structure  from  x-ray  dif- 
fraction measurements,  as  is  commonly  the  case,  has  required  sum- 
mation of  Fourier  series  of  the  general  form: 

®xyz  k 5 2«i(hx+ky+iz)  (l) 

where  the  sum  is  real  and  the  coefficients,  A,  h,  k,  are,  in 
general,  complex.  This  sum  may  be  written  in  the  form: 

zl^hke!  [2«(hx+ky+iz)  + . (2) 

In  order  to  obtain  sufficient  resolution  in  complicated  cases,  such 
a sum  may  require  evaluation  at  each  of  6o  values  of  each  of  the 
coordinates  in  the  range  0 - n.  Such  sums  are  necessarily  too  labo- 
rious unless  performed  by  punched-card  methods  or  by  suitable  pro- 
grammlng  on  high-speed  electronic  computers'  ' . (However,  the  very 

elegant  XRAC  cooputer  of  Pepinsky^^)  is  specially  designed  to  confute 
2-dimensional  sections  for  such  series.) 

In  the  earlier  stages  of  such  structure  investigations  simpler 
forms  of  the  series  are  often  useful.  These  may  often  be  1-  or  2- 
dimensional  and  centrosymmetric : 

®2  • fe  k ^(0  ^3) 
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and 


h *"h  *^00  2n  hx,  (4) 

and  the  coefficients  are  now  real.  For  the  purpose  of  computation  (3) 
may  be  expanded: 

cos  2n  kxj  cos  2n  hy 

-^|;Z  A^q  sin  2n  ltx|  sin  2n  hy.  (5) 


Evaluations  of  (4)  and  (^)  are  normally  carried  out  using  a manually 
operated  calculating  machine  and  Beevers-Llpson^**^  (or  Patter son- 
Tunell^^^)  strips  which  carry  values  A 

2rt 


rslnj 
(^cos ) 


2nhx  at  Intervals  of 


It  la  as  an  alternative  to  this  procedure  that  attention  was  turned 
to  the  possibility  of  devising  a simple  computer  to  perform  these 
summations.  To  be  useful,  such  a computer  must  have  great  ease  of 
operation,  require  little  or  no  maintenance,  and  furnish  results  of 
sufficient  accuracy  (l^  or  better). 

Numerous  approaches  have  been  made  to  this  problem  in  the  past. 
Huggins^^^  designed  a series  of  masks  on  photographic  film  with  proper 
distribution  of  density  to  permit  the  summation  of  (3)  directly. 

Woolf son^^)  has  replaced  the  photographic  film  with  masks  made  from 
sets  of  holes  drilled  In  brass  sheet.  Neither  of  these  methods  permits 
summing  over  indices  h and  k greater  than  about  10.  With  Huggins' 
method  it  is  difficult  to  get  very  good  accuracy. 

A novel  analog  machine  for  computing  a sum  of  the  type  (3)  has 
been  described  by  Shimizu,  Elsey  and  McLachlan^^\  This  machine  has 
been  built  with  eight  terms.  It  would  appear  to  be  cumbersome  and 
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quite  expensive  If  extended  to  many  more  terms.  However,  It  has  the 
great  advantage  of  Inherently  rapid  operation. 

Among  the  one -dimensional  machines  designed  previously  are  a 
number  of  different  types,  of  which  several  are  noteworthy.  The 
historical  prototype  and  well-known  harmonic  synthesizer  of  Michelson^^^ 
la  not  well  adapted  to  rapid  numerical  calculation.  More  recently 
designed  machines  involve  the  use  of  harmonically  tapped  potentlom- 
eters^^^^  and  mechMical  generators  based  on  the  principle  of  the 
scotch-yoke  crosshead^^^  which  was  also  incorporated  in  some  models 
of  Lord  Kelvin's  *T?ide  Predictors”. 

These  designs  all  have  certain  drawbacks.  Most  of  them  are  quite 
expensive.  Probably  the  most  attractive  one  from  the  standpoint  of 
cost  and  utility  is  the  machine  using  the  harmonic  potentiometer  but 
this  computer  is  reported  to  give  a gain  in  time  of  only  about  U;1 
over  the  use  of  strip-adding  machine  methods. 

Consideration  of  the  design  of  a computer  for  suoaiing  series 
of  the  type. 


y fa 


Bn  cos  2nnx 


(6) 


1 n ^bjj  sin  2jtnxJ^ 
requires  that  a decision  be  made  concerning  certain  main  features. 
Numerical  methods,  while  well  adapted  toflnal  calculations,  do  not 
appear  promising  from  a speed  or  cost  standpoint.  Analog  computers 
by  natiire  appear  better  suited  to  the  problem.  It  was  decided  to 
survey  the  methods  for  generating  and  adding  harmonic  functions. 
Formally,  such  a cooputer  would  have  the  components  shown  in  Fig.  1 


3 


FIGURE  1 

A]^  - - - - Aq,  Coeine  Generators,  cos  n 0. 
®1  " ■ " ■ Generators,  ein  n 0. 

” Sq,  Cosine  Amplitude  Controls. 

- - - - bjj.  Sine  Amplitude  Controls. 
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Briefly,  the  requirements  to  be  met  are  as  follows: 

(1)  The  cosine  and  sine  sections  should  have 
at  least  20  terms  each. 

(2)  The  generators  Ajj  (B^)  must  be  free  of 
overtones.  (Error  less  than  1^). 

(3)  Amplitudes  of  generators  must  be  readily 
controlled . 

(U)  A simple  method  of  addition  and  sub- 
traction must  be  available  so  that  the 
terms  can  be  added  or  subtracted  at  will. 

(5)  The  output  must  be  of  a type  easily  and 
quickly  read.  It  is  desirable  that  auto- 
matic recording  be  achieved  simply. 

(6)  Cost  should  be  low,  otherwise  punched- 
card  methods  compete. 

(7)  Over-all  simplicity  and  speed  of  oper- 
ation must  be  such  that  the  method  is 
a marked  improvement  over  the  strip- 
calculation  methods. 

Consideration  shows  that  the  generators  may  be  of  either  of  two 
classes,  namely,  (l)  those  for  which  the  amplitude  depends  on  the  time 
or  (2)  those  in  which  the  amplitude  is  a function  of  a positional 
coordinate  only.  The  first  class  is  exemplified  by  the  various  oscil- 
lators, electronic  and  otherwise,  in  which  the  amplitude  is  given  an 
expression  of  the  form 


A^  • A^  sin  k (7) 

For  the  purpose  in  hand,  such  generators  are  not  best  suited  since 
amplitudes  must  be  saddled  at  properly  chosen  times.  They  have  marked 
advantages  when  elaborated  to  the  extent  Pepinsky  has  gone  in  the  XRAC. 

As  a result,  we  have  confined  our  attention  to  generators  of  the 
second  class  of  which  the  "scotch-yoke"  and  sine-tapped  potentiometers 
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are  examples.  Such  generators  develop  amplitudes  which  are  functions 
of  a position  coordinate 

® (8) 

This  makes  for  marked  simplicity  in  withdrawing  a result  from  the 
machine . 

It  seemed  that  the  requirements  listed  might  be  met  best  by 
making  application  of  the  law  of  Malus  for  the  transmission  of  light 
through  a pair  of  polarizers: 

Itr  » Iq  cos^  e/2  » ^ (1+cos  9)  (9) 

Such  a generator  would  consist  of  a small  incandescent  filament 
light  source,  a rotatable  polarizer,  a fixed  polarizer  and  a photo- 
cell pick-up.  This  scheme  apparently  has  not  been  used  previously 
and  is  now  feasible  because  of  the  present  low  cost  of  polaroid 
material  and  of  commercial  phototubes.  The  maximum  amplitude  of  the 
term  is  controllable  by  adjustment  of  the  heating  current  for  the 
light  source.  Linearity  of  response  is  achieved  by  use  of  a vacuum 
type  photocell. 

A test  with  Type  "J"  Polaroid  in  conjunction  with  a 922  (S-1 
surface)  cell  indicated  that  equation  (9)  was  obeyed  within  ± 0.1^ 
except  for  an  additive  constant  attributable  to  a light  leak  in  the 
red.  On  recommendation,  Type  "H'*  Polaroid  was  selected  since  its 
predominant  light  leak  is  in  the  blue.  Transmission  measurements 
with  the  polarolds  in  crossed  position  and  a S-1  cell  still  showed 
an  appreciable  leak  in  the  red.  However,  use  of  a 929  (S-4  surface) 
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cell  reduced  the  response  to  light  leak  to  a negligible  value.* 

As  a result  of  these  testa  a design  for  a cooq;>lete  analog 
harmonic  synthesizer  was  worked  out.  Balancing  economy  against 
sufficient  flexibility  and  accuracy^  the  number  of  sine  and  cosine 
terms  was  set  at  twenty  each.  Figure  2 shows  the  generator  design. 

*nie  synchronized  drive  for  the  two  groups  of  generators  is 
similar  to  that  employed  by  Mlchelson  and  is  furnished  by  a cone  of 
spur  gears  each  having  lOn  teeth,  the  cone-shaft  being  driven  by  a 
60:1  worm  reduction  gear.  Bach  generator  carries  a 20-tooth  spur  gear 
to  engage  with  the  corresponding  cone  gear.  Thus  rotation  of  the  cone 
shaft  through  an  angle  of  9 rotates  a generator  through  an  angle  of 

n(^).  (Rotation  of  the  worm  shaft  by  one  revolution  is  equivalent 
to  6®  in  O).  Figure  3 shows  the  assembly. 

Figure  4 is  a photograph  of  the  assembly.  A pair  of  sine  and 
cosine  generators  are  on  opposite  sides  of  the  corresponding  cone  gear, 
supported  by  a cross  member  through  the  center  of  which  passes  the 
cone  shaft.  The  shaft  hole  in  each  case  is  bushed  with  a needle 
bearing.  In  the  lower  part  of  the  assembly  each  succeeding  cross 
member  is  advanced  through  an  angle  of  n/4  to  provide  clearance  for 
the  generators,  while  in  the  upper  portion  the  angle  of  advance  is 
necessarily  Increased  to  n/2. 

The  wiring  circuit  is  shown  in  Figure  Each  gang  of  light 
sources  is  powered  from  one -half  the  secondary  of  a six-volt,  center- 
tapped  transformer.  In  series  with  each  lamp  bulb  is  a 10-ohm  vlre- 
wound  rheostat  permitting  selection  of  (maximum)  amplitude. 

*We  are  indebted  to  Dr.  Donald  Tuoml  for  making  the  transmission 
measurements  and  suggesting  the  use  of  the  S-4  surface. 
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SECTION  X-X 


Ftg.  Z. 


T 


6 


9 


11 


t 


Addition  and  subtraction  are  siotply  effected  by  adding  the  photo- 
currents  in  parallel  or  antiparallel,  respectively.  For  this  purpose 
the  cosine  cells  (and  similarly,  the  sines)  have  as  voltage  supply  two 
135 -volt  B-batterles,  one  for  addition  and  the  other  for  subtraction. 
Each  pair  of  leads  from  a photocell  terminates  at  a 2-pole,  3-position 
switch.  The  three  iwsitions  correspond  to  + connection,  connection  to 
auxiliary  output  and  - connection. 

The  sum  actually  produced  is 

Sl‘  ®n  ®n  ®c  * f ^ f \ “ ®S 

in  which  the  coefficients  a^  and  b^  may  be  + or  - and  for  the  sine  terms 
®S  ■ ®C  phase  angle,  is  introduced  by  initial  adjust- 

ment of  the  fixed  polaroids  in  the  sine  generators.)  The  sum  required 
is 

Si  ■ a_  a cos  n 9 b_  sin  n 0 (6') 

y n n n n 

The  sum  is  obtained  from  thus 

Si  * Si'  - ir  \)  + a^,.  (11) 

The  removal  of  the  second  term  of  (ll)  is  accomplished  by  addition  (or 
subtraction)  of  a parallel  current  furnished  by  a 1.5-volt  dry  cell 
and  controlled  by  a separate  rheostat.  In  a similar  manner  a^  is  intro- 
duced. Arrangement  is  made  for  the  addition  or  subtraction  of  the  sine 
sum  and  the  cosine  sum. 

A view  of  the  control  panel  is  exhibited  in  Figure  6.  Shown  in  the 
upper  portion  of  the  panel  are  the  3-position  photocell  switches,  each 
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with  Its  anq;>lltude  control  knob.  The  set  of  20  cosine  terms  is 
Immedlatel}'  beloW  the  set  of  sines.  Belov  on  the  left  are  two  output 
Jacks.  The  upper  (jj^)  is  for  auxiliary  output  and  the  lower  (J2)  is 
a mlcroamtoeter  on  which  the  total  sum  is  read. 

Just  below  the  Jacks  is  a 4-position  switch  which  connects  to  the 
main  output  (l)^  cosines,  (2))_  sines,  (3)^  cosines  +)_  sines,  (4) 
y cosines  sines.  At  the  right  of  the  lower  panel  are  the  power 
switches,  and  switches  and  controls  for  introducing  and  removing  constant 
terms. 

In  practice  each  photocell  is  connected  in  turn  to  to  which 
is  connected  a sensitive  microammeter . The  control  for  the  term  in 
question  is  then  adjusted  to  produce  the  proper  amplitude.  This  is 
done  with  the  cone-shaft  worm  wheel  set  at  0 as  shown  on  the  revolution 
counter.  In  this  position  the  output  is  2 for  a cosine  term  and  bjj 
for  a sine  term  (by  equation  (lO)).  As  a result  twice  the  coefficient 
must  be  introduced  for  cosine  terms  while  the  coefficient  itself  is 
used  for  sine  terms. 

After  adjustment  of  these  coefficients,  the  terms  are  connected 
into  the  summing  circuit  + or  - depending  on  the  sign  of  a^^,  b^.  Now, 
with  the  4-position  switch  set  in  Position  1 the  cosines  is  deter- 
mined (i.e.  ^ 2 Sjj)  and  half  of  its  value  is  removed  by  adjustment  of 
the  proper  rheostat  at  the  right.  The  sines  are  treated  similarly, 
using  Position  2 of  the  4-position  switch,  except  that  the  output  is 
reduced  to  zero  by  the  proper  rheostat  adjustment  (at  the  right). 
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The  total  sum  can  now  be  read  on  the  main  output  in  Position  3 
or  4.  Each  half -revolution  of  the  worm  drive  shaft  corresponds  to 
change  in  0 of  3° • 

While  it  has  not  been  so  used,  a recording  microammeter  could  be 
connected  conveniently  to  the  output.  With  a synchronous  drive 
applied  to  the  worm  shaft,  each  sum  over  h at  constant  k of  (3)  would 
then  be  automatically  computed  as  a continuous  function  of  x.  The  sums 
over  k at  selected  x's  would  then  be  recorded  as  continuous  functions 
of  y. 

We  are  indebted  to  the  University  Committee  for  Allocation  of 
Research  Foundation  Grants  for  aid  covering  most  of  the  shop  work  of 
assembly.  We  wish  to  express  our  appreciation  to  Mr.  Gordon  Laverack 
and  Mr.  Audrey  Stennett  for  their  skillful  workmanship. 

Numerous  additional  applications  of  this  machine  are  possible. 

One,  of  interest  to  the  present  work,  is  in  connection  with  the  de- 
termination of  the  roots  of  a polynomial.  Since  it  is  possible  to 
take  out  sine  and  cosine  sums  separately  and  simultaneously  by  the 
introduction  of  an  extra  output  Jack  on  the  panel,  the  machine  will 
serve  as  an  'Isograph"  by  connection  of  these  outputs  to  an 
"x,y"  recorder. 

As  is  well  known  from  the  theory  of  functions  of  complex  variables, 

if  are  the  roots  of  a polynomial,  p(z)  * A_  for  p(z)  = 0, 

^ n»o  « 

then  the  number  of  roots  whose  magnitudes  are  less  than  an  arbitrarily 
chosen  number,  p,  can  be  determined  from  an  x,y  plot  of 

P-  - m 

X e p Aq  sin  ne  and  y a ^ cos  n0. 
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The  number  of  times  the  function  wraps  itself  around  the  point  x > 0, 
y • >Aq  is  then  the  number  of  roots  for  which  | z^j  is  not  greater  than 
p.  Iteration  then  will  lead  to  values  of  the  roots. 

This  problem  arises  in  the  direct  solution  of  the  structure  ampli- 
tude equation  for  coordinates  of  the  atoms  in  the  unit  cell,  for 
instance  as  outline  by  Avraml^^^).  For  an  especially  simple  (but 
frequently  met)  centro-symmetric  case,  the  amplitude  equation  becomes: 

^hOO  . 

^ * I cos  h . 

* atom  i*l 

It  is  here  assumed  all  atoms  have  identical  scattering  power,  f,  there 

“i 

are  N atoms  per  cell  and  the  parameters  ^ signs  as  well 

as  experimental  magnitudes  of  the  quantities  F^qq  are  presumed  known 
for  a sequence  of  N values.  These  may  be  readily  written  as  linear 
functions  of)  cos“  a^,)  cos  etc.  and  the  quantities  cos 

are  roots  of  the  polynomial 

2 3 N 

®o  ®1*  ®2*  ®3*  ■ — = 0. 

The  coefficients  a^  eu:e  easily  obtained  from  the  quantities 
Sjj  « L CO*  and  therefore  from  the  amplitudes,  Fj^qq. 
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